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Background 


The purpose of this research effort was to see if displacement or mixed 
version hp finite elements in time could be used to efficiently solve the 
rotor trim problem. At the end of the first year of effort it appeared the 
funding would end. Therefore, w f c prepared a "Final Report" that 
summarized the first year’s effort. The effort included all ol the linear 
results for flap and flap-lag with both displacement and mixed elements. 


Second-Year Effort 


However, it turned out that second-year funding did become available. 
Thus, the "Final Report" became a Semi-Annual Report. Therefore, the 
report given here summarizes the second year’s effort which is the solution 
of the nonlinear trim problem. As it turns out, this was also successful. 
The conclusion of our work, therefore, states that finite elements in time 
arc a good way to solve rotor trim problems for a moderate number ol 
degrees of freedom. The technical content of this second-year effort 
follows. 
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INTRODUCTION 


Finite Element m Tune has been proven to lie a powerful alternativ<’ strategy 
for rotor eralt trim problem. Recently. Finite Element Method in lime has been 
well developed in various versions like time-marchinc. framework. Calerkin framework. 
R a vlciuh- Kit z framework and mixed formulation. Hr. Peters and Mr. Hon applied 
this method on the rotorcraft trim problem to obtain linearized solutions and obtained 
a very good result. This project is aimed to expand the application in rotorcraft trim 
problem from linearized solution to nonlinear solution. 

The rotorcraft trim problem is to find a period solution for period-coefficient, differ- 
ential equations subject to side constraints that certain force and momentum balance 
equations have to be zero. There are certain free ( or trim ) parameters that must 
be chosen to meet these side constraints. A lot of successful works have been done 
by Dr. Peters and Mr. Hou in this area. They linearized the flap-lag equation and 
compared the results of Fourier series analysis, displace method and mixed method 
for accuracv and efficiency. Finite Element in Time presented a strong potential to 
outperform the conventional Fourier series analysis. Therefore, it is interesting to find 
out if Finite Element in Time can further provide a good flap-lag nonlinear solution. 
In this project, we use the displacement method for the hovering flight. 
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NONLINEAR EQUATION 
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The model which we use is the same as in Ref [1] [2]. The schematical 
model is shown in Fig 1. The reverse flow effect is neglected. We also 
assume that the drag coefficient is small with respect to the lift curve 

slope a. 

The euations of motion of this model are [2]: 
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For hovering flight 
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When 0 = 0 or /? = 0, there is no elastic coupling between flap 
and inplane. Then 

P= \ + co 2 p = p (o-o)^ Z - 0 
Expanding 0 as 

0= 0+ e p ( p- p pc ) + 0 ? C = ( 0 - Bp Ppc ) + + 0c C 

0 = 0o + 0y s/n f' + 0 C cos f* 

By small angle assumption 
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Substituting all these terms into equation (1), we obtain the 
nonlinear equations for flap-lag forced response 
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The constraint equations for trim condition are 
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When thrust is specified, equation (3) can be replaced by 
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For hovering, 


CK ' )= [ 


niSPl.ArF.MF.NT METHOD [1] 


For displacement method, we begin with the Hamilton's law of 
varying action 
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With a weak constraint of momentum, the hamilton's weak 
principle becomes 
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Where the trailing term 

) 

£ &?,/>, i;; 

1=1 

can be thought of as the virtual action 8A entering ( or leaving ) 
the system at time *i and U 

For a typical spring-damp-mass oscillator with unit mass 

L-T - V = J-.t 2 - \k x 1 
2 2 

SWnc = (F-CX)5X 

Then equation (8) becomes 

l ( x + cx+kx-F)5xdt - [(P_px(tf))Sx - ( P /-x(f,))5x]=0 (9) 

J ti 

) For this equation to be valid for any 8x , we must have 

x+cx+kx=F 

Pf = x(tf) 

Pi=x(ti) 

for periodic problem, ti= 0, tf=T, x(0)=x(T) is a strong condition 
and j s a condition enforced through Po = Pt ■ By 

introducing Lagrange multiplier, X , equation (8) becomes 

j ( xSx -Cx8x-kx5x + F8x)dt - Xp[ 8x(T) - &(0) ] = 0 (10) 

where Xp=P q=Pt 
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For the test and trial functions, there are certain geometric 
boundary condition which must be satisfied! 1 ] . 
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STRATEGY 


To solve flap-lag problem by displacement method, we use 
Galerkin's scheme, first, we choose the integrals ot legendre 
polynomials P/x) as shape functions for p, £, 8p and 5C, [ 4 ] 


£=X <t>j( i ) ) Qj 
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A. One Element Quasi-Trim Case 

For this case, the collective and cyclic pitch angles are given, 
$ = 0 O + Q^in t + G c cos t f so that Pc>Ps~® to minimize rotor hub 
moment. Multiplying equation (2) by (pi (i=l ,2 ,n) on both 
sides and integrating over [0,2 tt ] 

j [P <pi + C\ ip (pi + C 12 C <Pi + KnP <Pi K12 C 0 il d t = J F\<pidt 
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Set 


[COS 2 /^ (pi + C21/3 <pi + C22C + ^2lP 0 » + ^22 C 0 j] ^ * “ I ^ 20 i d 1 


P(2n) = l3(0) = Pp = Xp 


£( 2k )cos 2 fi ( 2tt ) = £( 0 )cos 2 (3 ( 0 ) = 


change variables 


r = tt( t+ 1 ) 


-1 <T< 1 


0=4— 1 4-= 1 o' 

dt K dx K 


After integrating by parts, we get 2n nonlinear equations of q\. 
Pi, Xp and Af (i=l, 2,--,n). The other two equations which are 
required to close the system come from the periodic constraint 

P(0) = P(T) C(0 ) = Z(T) 

n n 

=> £[0,(-l)-0/(l)]<7«=O 

i=l «=1 

The final system of nonlinear equations is 
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The matrix [A] includes /?, £, /T and £\ which are functions of 
unknown variables qi and pi. Therefore a iteration process is 
required. The first guess for qi and pi is not difficult. In fact, 
the stability of this system of nonlinear equations is 
tremendous. Even if the initial values of qi and pi. are set to be 
equal to 100 (although there is no physical meaning), they will 
converge to the correct result within a few steps. 


B. One Element Trim case 

In this case 9 0 9 c and 9 S are unknowns to be found such that (5 C 
=p s =0. So we have 2n+5 unknowns. Notice that in equation (2), 
the 9o,9 c and 9 S are included in {F}. We need to extract them 

from {F} and move them to the left hand side. The equations 
then become 

[m(A») 1 |?|+[c(j9,f.ftf.e)] J?j +[k(/},c,»)] j^j 

+ [t (9)1 So + sin t (T(0)) 9 S + cos t lr(0))0 c ={Pr} 

Where 
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\ Following the same scheme, we will obtain 2n+2 equations. The 

trim condition, provides the extra 3 equations to close the 
system. 

C. Multi Element Trim Case 

We divide [0,2 ti] into M elements, each has a length of 2rc/M, 
Within the m th element let 

p (m) = X <*>/ 77 ) q ) m) c (m) = X ( n )p) m) • ■> M 

j= 1 /=! 

By continuity condition 

By periodic condition 

Also we set the nodal momenta to be 

• (m) * (m-1) ( m \ 

P \t m ) = p (tm) = 4 ) 

m = 2,3,- • ., M 

C (,t m )COS 2 p K \t m )= C ( t m )C0S 2 p (t m ) = A£ 

P { \) = p (M \tM + l) = ^ 

£ \ti)C0S 2 p {l \t\) = ^ Vm+1 2 p (W) (tW+l ) = ^ * 

Then the trailing terms in equation (12) become 
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The trim constraint for multi elements with (N-l) th order 
polynomial are [1] 
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a cr 


The final system of nonlinear equations includes 2M(N+l) + 3 

, T , ^ m) , b< ot) , A« m) , A r m) , 00, 0c and 0 5 . 

unknowns. They are H t ,y i ' P s 


RESULT AND CONCLUSION 

The results are very close to those obtained by other methods. 
This varifies the validity of this method. As far as the flap lag 
problem is concerned, there is no restriction on the choice of 

1 0 



initial values. They always converge to the final results within 
a few iteration steps. Also notice that the results are almost 
independent of the number of elements and the order of 
Legendre polynomials. Therefore we can conclude that, in 
solving this kind of nonlinear problem, the displacement 
method has a very good stability. Since the flap lag problem in 
hovering flight is a very much simplified case compared with 
the flap lag problem in forward flight, additional work is 
required to investigate the performance of displacement 
method in forward flight condition. It is also a interesting topic 
to see how the mixed method behaves in such cases. 


no.Elements/pol.order 

p/ioo 

C/1000 

00 

e s /ioooooo 

e c /ioooooo 

1 

4 

9.6825644 

-6.3044904 

0.2970519 

1.9670248 

-1.7705892 

1 

6 

9.6825607 

-6.3045006 

0.2970519 

1.6492519 

-1.1021315 

1 

8 

9.6825600 

-6.3044997 

0.2970519 

1.6540506 

-1.0808242 

2 

4 

9.6825629 

-6.3045016 

0.2970519 

2.5363567 

-1.4900319 

2 

6 

9.6825629 

-6.3045030 

0.2970519 

1.1276192 

-9.1521546 

2 

8 

9.6825622 

-6.3045025 

0.2970519 

9.1574357 

-1.0367157 

5 

4 

9.6825719 

-6.3045076 

0.2970521 

1.0911533 

-9.4596531 

8 

4 

9.6825704 

-6.3045053 

0.2970520 

-6.8381580 

2.0519121 

8 

8 

9.6825697 

-6.3045057 

0.2970520 

-1.6251148 

-3.1796359 
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